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Level density of odd-A nuclei at saddle point 


Wei Zhang,!:* Wei Gao,! Gui-Tao Zhang,! and Zhi-Yuan Li! 
! School of Physics and Microelectronics, Zhengzhou University, Zhengzhou 450001, China 


Based on the covariant density functional theory, by employing the core—quasiparticle coupling (CQC) model, 
the nuclear level density of odd-A nuclei at the saddle point is achieved. The total level density is calculated 
via convolution of the intrinsic level density and the collective level density. The intrinsic level densities are 
obtained in the finite-temperature covariant density functional theory, which takes into account the nuclear 
deformation and pairing self-consistently. For saddle points on the free energy surface in the (2, y) plane, the 
entropy and the associated intrinsic level density are compared with those of the global minima. By introducing 
a quasiparticle to the two neighboring even-even core nuclei, whose properties are determined by the five- 
dimensional collective Hamiltonian model, the collective levels of the odd-A nuclei are obtained via the CQC 
model. The total level densities of the 7°4~?4°U agree well with the available experimental data and Hilaire’s 
result. Furthermore, the ratio of the total level densities at the saddle points to those at the global minima and 
the ratio of the total level densities to the intrinsic level densities are discussed separately. 
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I. INTRODUCTION 


As a fundamental nuclear property, nuclear level density 
(NLD) plays a crucial role in many applications, such as the 
calculation of reaction cross sections with nucleosynthesis, 
the nuclear reaction calculation program TALYS [1], and the 
Hauser—Feshbach model for compound nucleus alculation- 
s [2, 3]. Owing to the complexity of nucleon interactions and 
the fact that the level density increases exponentially with an 
increase in the excitation energy, the accurate calculation of 
the NLD has long been a theoretical challenge. 

Many methods for estimating NLD have been develope- 
d. The most common method is the Bethe formula based on 
the zero-order approximation of the partition function of the 
Fermi gas model [4, 5]. In attempting to reproduce the ex- 
perimental data, various phenomenological modifications to 
Bethe’s original analytical formulation have been proposed— 
particularly to account for shell, pairing, and deformation 
effects—which led first to the constant-temperature formu- 
lation, then to the shifted Fermi gas model, and later to the 
popular back-shifted Fermi gas model [2, 6, 7]. 

There are many microscopic methods for calculating NLD, 
including the shell-model Monte Carlo method [8-10], the 
moments method derived from random matrix theory and 
statistical spectroscopy [11, 12], the stochastic estimation 
method [13], the Lanczos method using realistic nuclear 
Hamiltonians [14], the self-consistent mean-field approach 
based on the extended Thomas—Fermi approximation with 
Skyrme forces [15], and the exact pairing plus independent 
particle model at a finite temperature [16-19]. On the ba- 
sis of the Hartree-—Fock—Bogoliubov (HFB) model, S. Hilaire 
and S. Goriely developed a microscopic approach to describe 
NLD with great success [20-22]. Microscopic methods based 
on the self-consistent Hartree-Fock (HF) plus BCS model 
[23-25] have also been developed. 

Recently, on the basis of the relativistic Hartree— 
Bogoliubov model [26-29], J. Zhao et al. developed a method 
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for describing NLD [30]. In this model, the partition func- 
tion is determined using the same two-body interaction em- 
ployed in the HF plus BCS and HFB mean-field models [25], 
which includes shell, pairing, and deformation effects self- 
consistently. The total level densities are the convolution of 
the intrinsic level density and the collective level density. The 
intrinsic level density is obtained by an inverse Laplace trans- 
form of the partition function with the saddle-point approx- 
imation [31, 32]. Previously, the collective enhancement is 
considered via a phenomenological or semi-empirical multi- 
plicative factor for rotational and vibrational degrees of free- 
dom [24, 33-36] or more microscopically via a combinatorial 
method using single-particle level schemes obtained through 
HF plus BCS or HFB calculations [20, 22]. In Ref. [30], 
the collective enhancement is determined from the eigen- 
states of a corresponding collective Hamiltonian that consid- 
ers quadrupole or octupole degrees of freedom. Both the in- 
trinsic level density and the collective enhancement are de- 
termined by the same global energy density functional and 
pairing interaction. 

The success of the microscopic description of NLD in 
even—even nuclei prompts us to extend it to odd-A nuclei. 
The core—quasiparticle coupling (CQC) model introduces a 
quasiparticle to the neighboring even—even core nuclei within 
the same covariant density functionals and achieved progress 
for describing the quantum phase transition in odd-A nu- 
clei. It is based on the covariant density functional theory 
(CDFT), which has achieved considerable success in describ- 
ing ground-state properties of both spherical and deformed 
nuclei all over the nuclear chart [26, 27, 37—40]. Its successful 
applications include superheavy nuclei [41—45], pseudospin 
symmetry [40, 46, 47], single-particle resonances [48-50], 
hypernuclei [51-56], thermal shape transition [57-59], and 
shell correction [60—64]. 


In this study, the CQC model is applied to the calcula- 
tion of collective levels in even—odd uranium isotopes. For 
even-—even isotopes, the collective levels are attained via the 
five-dimensional collective Hamiltonian (SDCH) model [65]. 
Similar to Ref. [30], the intrinsic level density is obtained us- 
ing the finite-temperature CDFT [58, 59, 66]. Because the 


77 level density of the saddle point plays an important role in 
7s the compound nuclei reactions [67, 68], the level density of 
79 the saddle point psq and the level density of the global min- 
s ima Pmin are analyzed. The ratio of the level density at the 
s1 Saddle point to the level density at the global minimum is the 
s2 quantity of interest. 

ss The remainder of this paper is organized as follows. The 
theoretical framework is introduced in Sec. II. The results for 
es 734-2401) are presented in Sec. III. Sec. IV presents a short 
s summary. 


Il. THEORETICAL FRAMEWORK 


87 


ss With the assumption of decoupling between intrinsic and 
s9 collective degrees of freedom, the excitation energy of a nu- 
90 cleus can be written as E* = FE; + Ee, where Ee represents 
o1 the collective excitation energy. The total level density is ob- 
92 tained as [36] 


rs Prot (E*) = J 0: (Ei) pe(E" -Ejda (1) 


with the collective level density given as 


95 pe(E) = X8 (E — Ee) Te (Ec) . (2) 


æ For a collective state with the angular momentum Te, the de- 
generacy is Te (Ee) = 21e + 1. 

, ss The intrinsic level density p; can be obtained from the gi- 
ə ant partition function of two types of particles via the inverse 
Laplace transform and saddle point approximation [31, 32]: 
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P= AD (3) 


101 
where S represents the entropy, and D is the determinant of 
a 3 x 3 matrix that contains the second derivatives of the en- 
tropy with respect to the inverse temperature 8 = 1/(kpT) 
and u, = Bà (T = N, Z) at the saddle point. The intrin- 
sic excitation energy is calculated as £;(T) = E(T) — E(0), 
where £(T’) represents the binding energy of the nucleus at 
temperature T. The specific heat is defined by the relation 
Cy = OE; (T)/OT. 

According to the ideas presented in Ref. [32], the determi- 
nant D can be simplified to the following form: 
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where N and Z represent the numbers of neutrons and 
protons respectively, and A,(7 = N, Z) denotes the neu- 
tron (proton) Fermi surface. For convenience, the tempera- 
ture used is kgT (in units of MeV) and the entropy used is 
S/kp (dimensionless). 

Entropy is extracted in the finite-temperature covariant 
density functional theory. In the covariant density functional 
theory, the Dirac equation for single nucleons is 


[u i0" — V”) — (m + S)] Yk = 0, (5) 


122 where m represents the nucleon mass, and y(r) denotes the 
123 Dirac spinor field of a nucleon. The scalar S(r) and vector 
124 potential V” (r) are 

S(r) =asps + Bsps + Ysp3 + 8sAps, 
V” (r) =avjý + wv Gb) + ôv Ajy 


+ T30TV Jpy + Tarv Ajhy + eA", 


(6) 


12 respectively. The isoscalar density pg, isoscalar current j{,, 
127 and isovector current jee have the following forms: 


ps(r) =X pelr yrl) [vg (1 — fe) + Fe] 
k 


128 jy(r) = DD pelr)” ve(r) [uk (1— 2fk) + frl (7) 
k 


Fey (r) = >> balr)Py"be(r) [vg (1 — 2fe) + fe], 
k 


122 Where fp represents the thermal occupation probability of 
130 quasiparticle states, having the form fg = 1/ (1 + e¥r/kBT), 
131 Ek represents the quasiparticle energy for the single-particle 
i state k, and Ep = [(e, — A)? + A?]2. The BCS occupation 
133 probabilities v? and the related u? = 1 — v7 are obtained as 


134 follows: 
1 Ek — À 
2 
= — 1 — 
Uk 9 ( Ex ) , m 
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Uk = 2 E; 


136 A; is the pairing gap parameter, which satisfies the gap equa- 
137 tion at a finite temperature: 
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139 The particle number N(Z) is restricted by N(Z) 
2 > [v2 (1 — 2fk) + fk]. The entropy is computed using the 
k>0 


Ak (9) 


relation 


S = -kg X [fkn fe+(1— fr) In (1 fr)]- 


k 


(10) 


For even—even nuclei, the collective levels are obtained vi- 
a the five-dimensional collective Hamiltonian (SDCH) [65]. 
All the collective parameters, such as the inertia parameter- 
s and the collective potential, are extracted from the con- 
strained CDFT+BCS in the triaxial deformation space. 
143 For odd-A nuclei, the collective levels are calculated us- 
149 ing the CQC model [69], whose collective Hamiltonian is ex- 
150 pressed as 


143 
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where À denotes the Fermi surface, and £4^+*! and H4+! 
represent the single-particle energy and the collective ex- 
citation energy for the even—-even A + 1 core, respective- 
ly. T and A denote the mean and pairing fields associated 
with long-range quadruple—quadruple particle-hole interac- 
tions and short-range monopole particle—particle interactions 
between the odd nucleon and core, respectively. The I field 


can be expressed as 
p 2 7 F 
y d a F Z Hos 
+1 


p4+ 1 = 


ss ka A+1 
Qollu j 


x (en ajr 
(12) 


are the 


7 A+1 P A 
vrer (2 0)" malnove) 
reduced quadrupole matrix elements of the spherical hole 
(particle) and cores, respectively. The Fermi surface \ and 
coupling strength of the quadrupole field x are left as free 
parameters that are fit to data separately for positive- and 
negative-parity states. 

In this study, the 5DCH and CQC models are based on the 
CDFT calculation with the harmonic oscillator basis N¢ = 
16. 


Ill. RESULTS AND DISCUSSION 


The parameter sets of covariant density functional theory 
used in this study are PC-PK1 [70] and DD-LZ1 [71]. PC- 
PK1 is one of the most widely used point-coupling parameter 
sets, and DD-LZ]1 is a density-dependent parameter set that 
aims to alleviate the spurious shell closure. The pairing effect 
is considered by the separable pairing force [72]. The nuclei 
considered are even—-even 734—74°U and even—odd ?35-239U, 

In the first step, large-scale finite-temperature CDFT calcu- 
lations are performed for 2°4~24°U in the temperature range 
of 0-2 MeV in the (2,7) plane. Fig. 1 shows the free ener- 
gy surface evolution with respect to the temperature for ?35U. 
The deformations of the global minimum and the saddle point 
change slightly with an increase in the temperature; i.e., with 
the increase in the temperature, the global minimum deforma- 
tion ə decreases, while the saddle point deformation y slow- 
ly moves toward the prolate axis (y = 0°) and G2 remains n- 
early constant. The saddle point gradually becomes indistinct 
as the temperature approaches 1.6 MeV. The free energy sur- 
face evolution for other nuclei 234-236-2400 is similar to that 
for ?35U. For the parameter set DD-LZ1, the energy surface 
shapes are similar, while the fission barrier heights are larger. 
Moreover, the free energy surfaces of ?34-240U considering 
the axial octupole deformations are checked, and it is found 
that although the PESs of some nuclei are relatively soft in the 
octupole direction, there is no significant octupole deforma- 
tion for the global minimum and first saddle point. Therefore, 
in the following calculation, the deformation space is limited 
to the (32, y) space. 

It is shown that the intrinsic level density has an expo- 
nential relationship with the entropy in Eq. (3). To study 
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Fig. 1. (Color online) Free energy surfaces in the (82,y) plane at 
temperatures of (a) 0, (b) 0.4, (c) 0.8, (d) 1.2, (e) 1.6, and (£) 2.0 
MeV for 72°U obtained via finite-temperature CDFT calculations 
using the parameter set PC-PK1. The global minimum and saddle 
points are represented by squares and stars, respectively. The energy 
separation between contour lines is 0.5 MeV. 


the dependence of the entropy on the nuclear deformation- 
s, Fig. 2 shows the entropy surfaces at T=0.4 MeV and 0.8 
MeV in the (2, y) plane calculated using two parameter set- 
s PC-PK1 and DD-LZ1. The global minima and the saddle 
points displayed in free energy potential surfaces Fig.1 are 
also marked in Fig. 2. Comparing these two figures reveals 
that the entropy has a low value near the free energy global 
minimum, and it becomes high near the saddle point. There 
are numerous similarities between the entropy surface and the 
free energy surface. For the low temperature of T=0.4 MeV, 
there are sharp changes within a certain deformation range— 
particularly for DD-LZ1. For the high temperature of T=0.8 
MeV, the entropy surfaces share substantial common features. 
Furthermore, several derivatives composing the VIDI term, 
which appears as the denominator in Eq. (3), are extracted for 
all the deformation grid points, and the intrinsic level density 
pi is settled. The intrinsic level density surface on the loga- 
rithmic scale is analogous to the entropy surface here. It is 
omitted owing to space limitations. 

In the nuclear reactions, the level density of the saddle 
point is critical. Fig. 3 shows several properties of the saddle 
points with respect to the temperature for 2°4~4°U calculat- 
ed with the parameter set PC-PK1, i.e., the excitation energies 
E;, the entropy S, the specific heat Cy, the partial derivative 
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Fig. 2. (Color online) Entropy S for 5U at T = 0.4 MeV (upper 
row) and 0.8 MeV (lower row) obtained via CDFT calculations using 
the parameter sets PC-PK1 (left) and DD-LZ1 (right) in the (62,7) 
plane. 
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Fig. 3. (Color online) Excitation energy Æ; (a), entropy S (b), spe- 
cific heat Cy (c), 0S/OT (d), pairing energy gap A (e), and intrinsic 
level density pi (£) with respect to the temperature for saddle points 
of ?34—240U, The results were obtained via finite-temperature triax- 
ial CDFT calculations with the parameter set PC-PK1. 


25 of the entropy with respect to the temperature S/T, the 
226 pairing energy gap A, and the intrinsic level density p;. The 
2237 temperature range ends at 1.4 MeV, which corresponds to the 
22 excitation energy of approximately 40 MeV. For higher tem- 
22 peratures up to 2 MeV, the trends do not change. In Fig. 3(e), 
230 pairing phase transitions occur at a temperature of approxi- 
231 mately 0.6 MeV for all nuclei. When T = 0.4 MeV, more 
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Fig. 4. (Color online) Same as Fig.3 but for the global minimum. 


232 curve details about Cy and 0.$/0T in Fig. 3(c) and (d) can 
233 be observed. For Cy, the curve reaches a local maximum and 
234 then decreases slightly, and the segments T < 0.4 MeV and 
25 T > 0.6 MeV are essentially two straight lines. Because the 
specific heat is the partial derivative of the excitation ener- 
gy with respect to the temperature, it supports the fact that 
the excitation energy curves in Fig. 3(a) are actually two con- 
necting quadratic parabolas. This is consistent with the fact 
that in the Fermi-gas model, the intrinsic excitation energy in- 
creases quadratically with respect to the temperature E; x T? 
with slope changes around the pairing phase transition. For 
OS/OT, it reaches 60 MeV~", decreases to 40 MeV~!, and 
then remains constant for T > 0.6 MeV. Because of the di- 
rect relationship between 0S/0T and entropy S, the entropy 
should be quadratic for T < 0.4 MeV and linear for T > 0.6 
MeV. This confirms the classic relation S œ T but only for 
high temperatures. The logarithmic intrinsic level density pi 
in Fig. 3(f) has analogous temperature dependence to the en- 
tropy S in Fig. 3(b), representing In(p;) x S in Eq. (3). 

For comparison, the corresponding properties of the global 
minimum of ?34-240U are shown in Fig.4. In Fig.4 (e), the 
pairing phase transition occurs at T, = 0.6 ~ 0.7 MeV, while 
254 the proton gap is larger than the corresponding neutron gap at 
255 low temperatures. For the specific heat Cy in Fig.4 (c), the s- 
256 lopes before and after the phase transition are close, rendering 
257 the excitation energy curve in Fig.4 (a) a smooth parabola. At 
25s high temperatures, the partial derivative 0S/OT in Fig.4 (d) 
259 has nearly the same constant as that in Fig.3 (d). Because the 
260 intrinsic level density increases exponentially with respect to 
261 the entropy, this indicates that both the entropy S$ in Fig.4 
22 (b) and Fig.3 (b) and the logarithmic intrinsic level density in 
23 Fig.4 (f) and Fig.3 (£) have the same trends at a high excitation 
264 energy. 


236 


265 


26 


D 


26 


=] 


26 


© 


269 


27 


=] 


27 


272 


27. 


Q 


Í 274 
27 


a 


276 


27 


qq 


27: 


@ 


27! 


o 


281 


te] 


281 


band1 band2 band3 
—— —— —o- PC-PK1 


—a— —ẹ— —e— Exp. 
—>— —o— DD-LZ1 


a 


E, (MeV) 


Fig. 5. Calculated low-energy positive-parity bands (panels a and 
b), of 2°°U based on triaxial CDFT calculations with the parameter 
sets PC-PK1 (red) and DD-LZ1 (blue) together with the available 
experimental data (black) [73]. 


In addition, the result for the parameter set DD-LZ1 is ob- 
tained. The specific heat Cy varies more gently near the pair- 
ing phase transition, while other properties are roughly the 
same as those for PC-PK1. 

In the second step, the collective level densities of odd-A 
nuclei are calculated via the CQC model, while those of even— 
even nuclei are obtained via the SDCH model. The collective 
Hamiltonian of the odd-A nuclei is obtained by coupling the 
core parts of the two adjacent even—even nuclei and one par- 
ticle or hole in the spherical case. The two free parameters 
of the model, i.e., the Fermi surface ÀA and coupling strength 
x, are adjusted according to the experimental values of the 
low excitation spectrum. Details can be found in Ref. [69]. 
The CQC parameters for the odd-A ?35-239U nuclei corre- 
sponding to two density functional parameter sets PC-PK1 
and DD-LZ1 are presented in Table 1. 


TABLE 1. CQC parameters for 79°:?3?3°U based on triaxial CDFT 
calculations with the parameter sets PC-PK1 and DD-LZ1; the units 
of À and x are MeV and MeV/b, respectively. 


density functional nucleus parity A x parity A x 
35 + -5740 - -6911.5 
PC-PK1 237TU + -7140 - -6.7 11.5 
23 + -6.140 - -5.2 85 
2350Y 4 -7940 - -7.2 11.0 
DD-LZ1 37 + -5.07.0 - -6911.0 
23 + -6.040 - -8.1 8.0 

Take ?35U as an example. Its low excitation spectra ob- 
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tained using PC-PK1 and DD-LZ1 together with the experi- 
mental data from the NNDC [73] are shown in Fig. 5. The 
coupling strength and Fermi surface are finely tuned to re- 
flect the collective enhancement induced by negative parity 
according to the corresponding experimental data of the low 
excitation spectrum. The calculated levels exhibit good qual- 
itative agreement with the experimental results. 

In the final step, when performing the convolution using 
the Eq. (1) to obtain the total level density, it is found that 
the total level density pto depends on the number of collec- 
tive levels considered. This implies that a specific collective 
truncation parameter should be introduced and adjusted man- 
ually. To alleviate this problem, inspired by Ref.[36, 74], a 
factor exp(— E/T) is inserted into Eq. (1), and the total lev- 
el density is rewritten as 


E* — E; 
por(E*) = f exp(-= 


)pi (Ei) pc (E* — Ej) dE. 

(13) 
This factor simulates the coupling between the collective lev- 
els and intrinsic levels: higher collective levels correspond to 
a weaker contribution to the convolution; thus, a fixed trun- 
cation parameter is not needed. According to the highest ex- 
citation energy considered, tests reveal that an angular mo- 
mentum up to 18 A with the first 20 collective levels for each 
angular momentum is sufficient. 

In Fig. 6, the total level density and the intrinsic level den- 
sity at the global minima and saddle points are compared 
with available experimental data below 5 MeV [75], available 
calculations in the TALYS-1.95 package [1] directory den- 
sity/ground/hilaire or density/fission/hilaire/Max1 and direc- 
tory density/ground/goriely or density/fission/goriely/inner, 
or the Reference Input Parameter Library (RIPL) [20, 22]. 
Clearly, the intrinsic level densities deviate from the experi- 
mental data, indicating the necessity of collective degrees of 
freedom. For excitation energies up to 10 MeV, the total lev- 
el density calculated via the two parameter sets PC-PK1 and 
DD-LZ1 agrees well with the two sets of theoretical calcu- 
lations performed by S. Goriely and S. Hilaire in TALYS. 
However, for energies above 10 MeV, the total level densi- 
ty obtained via the relativistic density functional agrees well 
with Hilaire’s calculation, whereas Goriely’s calculation ex- 
hibits a sharp increase. The ratio of the total level density 
to Hilaire’s calculation result ptot/PHilaire is well within an 
order of magnitude. 

Furthermore, we compare the total level densities at the 
saddle points psq and those at the global minima pmin and 
plot their ratio in Fig. 7. For low excitation energies up to 5 
MeV, the majority of psq values are smaller than pmin, which 
may be related to pairing association [76]. The pronounced 
peaks at 0-5 MeV for 2°7~299U are actually caused by a small 
drop in the intrinsic level density at the global minimum (Fig. 
6). For intermediate energies in the range of 5 MeV < E* < 
30 MeV, psa increases more quickly than pin. In particular, 
for ?35U, this ratio increases linearly with respect to the exci- 
tation energy. When the atomic number increases from 7°4U 
to 2*°U, the curve of this ratio tends to shift downward. 

Finally, it is convenient to study the ratio of the total 
level density to the intrinsic level density K.on(E*) 
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Fig. 6. (Color online) Total level densities (solid) and intrinsic level 
density (dash-dotted) at the global minima (a-g) and saddle points 
(a'-g') of the 734~?4°U with respect to the excitation energy E*. 
The pink and blue lines correspond to calculations with the param- 
eter sets PC-PK1 and DD-LZ1, respectively. The experimental data 
(black squares) from Ref. [75] and calculations from Refs. [20, 22] 
are compared. 


Ptot(E*)/p;(E*) [30] as a collective enhancement factor and 
attribute it to the inclusion effect of collective levels. Fig. 8 
shows this collective enhancement factor Ko); with respect to 
the excitation energy E* for 2°4~4°U. Clearly, the curves can 
be divided into two groups: one group consists of low-lying 
concentrated curves for even-even nuclei 2°4~24°U, and the 
other group consists of scattered curves for three odd-A nu- 
clei ?35-239U, This implies different collective spectrum pat- 
terns for even—even nuclei and odd-A nuclei. For even—odd 
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Fig. 7. Ratio of the total level density at the saddle point to that at the 
global minima psa /Pmin in the logarithmic scale with respect to the 
excitation energy for 7°4~?°U based on triaxial CDFT calculations 
with the parameter set PC-PK1. The embedded subfigure is for 7?°U 
in the linear scale. 
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Fig. 8. Collective enhancement factor Keo with respect to the 
excitation energy for global minima of 7°4~?4°U based on triaxial 
CDFT calculations with the parameter sets PC-PK1 and DD-LZ1. 


nuclei, the factor Ko), obtained via the parameter set DD- 
LZ1 exceeds that for PC-PK1. The magnitude of Keo) is sim- 
ilar to that of the nuclei 94:9%98Mo, 106-108Pd, and 196:112Cq 
discussed in Fig. 5(a) of Ref. [30] at the low excitation en- 
ergy. The collective enhancement factor shown in this figure 
undoubtedly indicats the fact that the collective levels of odd- 
A nuclei are lower and denser than those of even-even nuclei. 
For level densities at saddle points, this factor Keo is almost 
indistinguishable from those of the global minimum. 


IV. SUMMARY 


In this study, according to the finite-temperature covariant 
density functional theory, the total level densities were ob- 
tained by convolving the intrinsic level densities and the col- 
lective levels achieved using the parameter sets PC-PK1 and 
DD-LZ1. For saddle points on the free energy surface in the 
(32,7) plane, the entropy and several derivatives composing 
the VID] term were extracted, and the intrinsic level density, 
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